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Q ■ Abstract 

We obtain the effective gravitational equations on the codimension-2 and codimension-1 branes in the 

in ■ 

cascading gravity model. We then apply our formulation to the cosmological case and obtain the effective 
Friedmann equations on the codimension-2 brane, which are generically given in terms of integro-differential 
equations. Adopting an approximation for which the thickness of the codimension-2 brane is much smaller 
Oh) than the Hubble horizon, we study the Minkowski and de Sitter codimension-2 brane solutions. Studying 

D 

^ (— I the cosmological solutions shows that the cascading model exhibits the features necessary for degravitation 

of the cosmological constant. 

> 1 1 Introduction 

On 
00 . 

, The cosmological constant problem is one of the most pressing conceptual problems in physics. This problem 

arises because the observed value of the vacuum energy is very small as compared to the values inferred from 
quantum field theory. Recently, a braneworld model which could provide a promising framework for addressing 
the cosmological constant problem has been developed. This model is the so-called cascading gravity model 
|TJ HI [31 H] , which could induce infrared modifications of gravity that can screen the effect of a cosmological 
constant. This idea, referred to as 'degravitation' [I], could possibly provide a dynamical solution to the 
cosmological constant problem, since any large cosmological constant initially present would degravitate away 
over time. 

In the cascading gravity model which is a generalization of the Dvali-Gabadadze-Porrati (DGP) model [5] 
to higher dimensions, one constructs a sequence of branes with decreasing dimensions placed one on each other, 
where each brane action contains the induced gravity term. The gravitational force falls off faster in large 
distances in the cascading model than in the original five-dimensional (5D) DGP model [2j. In the simplest 
six-dimensional (6D) cascading model, our 4D Universe, codimension-2 brane, is placed on a codimension-1 
brane which is embedded into a 6D bulk spacetime. The 6D model contains two cross-over scales, r% := jp? 
and r4 := -t-ttj where Mq, M§ and Mi are gravitational energy scales in the bulk and on the codimension-1 and 
codimcnsion-2 branes, respectively. Assuming that <C r^, it is expected that the gravitational potential on 
the codimension-2 brane cascades from the 4D regime at short scales, to the 5D one at intermediate distances 
and finally to the 6D regime at large distances. This model addresses a problem in the 6D brane world models 
with the induced gravity. If there is no induced gravity term on the codimension-1 brane, the bulk graviton 
propagator diverges logarithmically near the position of the codimension-2 brane. Then the energy scale r^ 1 acts 
as an infrared cut-off for the propagator [3] so that it remains finite even at the position of the codimension- 
2 brane. The cascading model suffers a ghost instability if there is no tension on the codimcnsion-2 brane. 
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However, very interestingly, it has been shown that there is a critical tension above which the model becomes 
ghost free [H [5] ■ This model may also provide a mechanism of the degravitation [5J [5] which can support a 
very small expansion rate of our Universe even in the presence of a large cosmological constant, namely the 
codimension-2 brane tension. Hence, the degravitation could provide a way to resolve the cosmological constant 
problem. 

A purpose of this paper is to see whether in reality there could be some features of degravitation in the 
cascading gravity model. In order to establish if this model exhibits degravitation it is necessary to understand 
its cosmological evolution and obtain the effective Friedman equations. The cosmological behavior can be very 
different from the ordinary 4D cosmology and the 5D DGP model, although they should be recovered in a 
certain limit. Cosmology in the cascading gravity model has been studied in the context of the 5D theory in 
[7], which is composed of gravity coupled to a scalar field [§], originated from the bending of the codimension-1 
brane in the 6D bulk. This theory is a nonlinear extension of the weak gravity limit of the full 6D model, but it 
is not the unique extension. Although the 5D theory may possess several similarities to the full 6D theory, the 
final confirmation should be made in the context of the original 6D theory. A regular solution involving a flat 
codimension-2 brane with a tension has also been discussed in [S]. Some implications of the cascading gravity 
model to cosmological events have been discussed in, e.g., Ref. (TOl [TTj . 

In this paper, we present the covariant formulation of the nonlinear effective gravitational theory on the 
boundaries, codimension-1 and codimension-2 branes. We then apply them to find the cosmological solutions, 
in particular the de Sitter solutions which may possess features of degravitation. In the gravitational equations 
on the codimension-2 brane, the bulk contributions are given in terms of the integration over the sixth direction 
where we take a finite thickness of the codimcnsion-2 brane into consideration. We study some cosmological 
solutions that show degravitation of the cosmological constant in the cascading model. In these solutions the 
effect of the effective vacuum energy into the Hubble expansion rate is suppressed by the ratio ^ in comparison 
with the naive expectation from the ordinary 4D cosmology. Moreover, we will discuss some self-accelerating 
solutions by applying the small thickness approximations where the codimension-2 brane's thickness is assumed 
to be much smaller than the size of the Hubble horizon. 

2 Gravitational equations in the cascading gravity model 

In the simplest realization of the cascading gravity our codimension-2 brane world is placed on a codimension-1 
brane embedded into a 6D bulk space-time. We start from the general ADM form of the 6D metric, with the 
bulk coordinate y playing the role of a time variable 



{6) ds 2 = N 2 dy 2 + g ab (dx a + N a dy)(dx b + N b dy), 



(1) 



1 The 6D metric is 




N 2 + N a N a TV, 
N a g a 
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Figure 1: The configuration of the codimension-f (plane) and codimcnsion-2 (thick line) branes in the six- 
dimensional bulk is shown. The induced metric on the codimension-1 brane is g a b and on the codimension-2 
brane is <7 M „. 



where g a b is the 5D induced metric on the codimcnsion-1 brane with a,b = z, fi (here Greek letters denote the 
4D space indices.). Having the bulk and brane coordinate systems, the full 6D action can be written as 



mat 
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S = My d 6 x^Tg^R + J d 6 x^Tg{^- {5) R + C 

+ J ( fxV^'g{^- {A] R + Cr t )8{y)8{z). (2) 

The codimension-1 brane is located at y — and the codimension-2 brane is placed on the codimcnsion-1 brane 
at y = z = (see Fig. 1). We are interested in derivation of the gravitational equations on the codimension-2 
brane. In the ADM formalism, the full 6D action for the cascading setup, including the appropriate Gibbons- 
Hawking boundary term in 6D part, can be rewritten as 



S 



+ JcfxV^Tg^^R + Cr^Siz), (3) 
where the extrinsic curvature K ao is given by 

Kab = ^{dygab ~ V a iV 6 - V b N a ), (4) 

and K = K a bg ab . To derive the gravitational equations we should take the variation of the action J3]) with 
respect to the variables N, N a and g ab which vanish at infinity but are nonzero on the boundary branes. Thus the 
boundary terms arising in the variations give the contributions to the gravitational theory on the hypersurface 
of y = (codimcnsion-1 brane) and hypersurface of y — z = (codimension-2 brane). In the following sections, 
we consider the variation of the 6D, 5D and 4D parts of the action and derive the effective gravitational theory 
on the codimcnsion-1 and codimension-2 brane, separately. 
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2.1 Variation of the 6D part of the action 

Taking the variation of the bulk terms with respect to N, N a and g ab gives the bulk Einstein equations 

^G AB =0, A,B = y,z,n. (5) 

The boundary contributions in the variation of the 6D part of the action give the contributions localized to 
the hypcrsurface y = 0. Since there is no term coming from 5d y N and Sd y N a the only terms localized to the 
codimension-1 brane come from the term proportional to 6d y g a b- Since 5g a b does not vanish on the brane, the 
boundary contributions to the variations of the 6D part of the action read 

- M Aj/ J dzd A x^)~g (Kg ab - K ab )Sg ab . (6) 



Here we have used y—^g = Ny—^g and A y is the discontinuity in a given quantity A over the codimension-1 
brane namely, A y A = 2^41^=0+. From now on we omit the subscript | 0+ - 

2.2 Variation of the 5D part of the action 

Now we take the variation of the 5D part of the action (J3j) with respect to g ab which gives us the contribution of 
the 5D part on the codimension-1 brane. Since there is a codimension-2 brane embedded on the codimension-1 
brane, there are also contributions localized to the hypcrsurface y = z = 0. To obtain the 5D contributions 
localized to the codimcnsion-2 brane, we go to the ADM form of the 5D metric g a b adapted to the codimension-2 
brane where z plays the role of a time variable^) 

( 5 W = g ab dx a dx b = N 2 dz 2 + g^(dx^ + N^dz){dx v +N v dz), (7) 

where g^ v is the 4D induced metric on codimcnsion-2 brane. In the ADM formalism in the 5D spacctime, the 
5D part of the action ((2]) can be rewritten as 

dW^(^ (5) #+-cr')<%) = I 'dW^( (4) i?+/c 2 -/c,,/cn%) 



where 



cfxV^gCf^iy), (8) 



1C ' tu= 2A7( 5 ^v-V M AC- V„Ag, (9) 
and JC = g^K.^ . Varing the 5D part of the action with respect to g a b leads to 

i J dzd i x^/^)~g{Ml ^G ab - S ab )6g ab , (10) 

which is localized to the hypcrsurface y — 0. Here S ab is the energy momentum tensor for the matter on the 
codimension-1 brane. The contributions localized to the codimension-2 brane come from the variation bd z g^ v . 
Note that there is no term which is proportional to 8d z N^ and 5d z M in the variation of the 5D part of the 



2 The 5D metric is 

9ab — [ .. 
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action. Thus the only boundary term localized to the codimension-2 branc coming from the variations of the 
5D part of the action is 

]\/[3 r 

- -^A 2 / d+x^/^TgiW - JCnSg^- (11) 

where we have used y— = My—^g. The discontinuity in a given quantity A across the hypersurface of 
z = is given by A Z A = 2^4| 2= o+. From now on we omit the subscript | 2 =o+- 

2.3 Variation of the 4D part of the action 

Finally, varying the 4D part of the action with respect to g ab leads to 

l - J d 4 x^r g (Mi (4) Gfit, - T^)5^ b 5g ab , (12) 
where T^ v is the energy momentum tensor for the matter on the codimcnsion-2 branc. 



2.4 Effective gravitational equations on the branes 



Collecting the boundary contributions in the variations of all parts of the action ([2]) derived in the previous 
sections, one can obtain the effective gravitational equations on the codimcnsion-1 and codimcnsion-2 branes. 
Combining (|6]) with (fT0|). and (fl2"j) . the boundary equations on the codimension-1 brane at y = read 



dzd^V^g (Ml^G ab - Sab) 



A, 



dzdW-^S M*(K ab - Kg ab ) 



(13) 



+ 



dzcfixV^Wg - Ml^G^ + + 2M 5 3 (/C„ 1/ - Kg^) 6£%6(z). 



Now we assume that the shift vector N a = n a (x)s(y)e(z). The function e(z) is a regulating function with the 
following properties: e(oo) = 1, e{—z) = — e{z) and e(z) x = 25 £ (z), where S e (z) is a regularization of the Dirac 
delta function. The function s(y) is the sign function and s(y)^ y = 25 (y). Using this ansatz one can see that 
the first term on the right hand side of equation (fl"3]l includes some terms localized to the codimension-2 brane. 
One can see this using the components of the extrinsic curvature K ab 

1 



K^u = — 

K zz 
K 



L Zfl 



2N 

1 
2N 

1 
2N 



d y g^ - d {fl N v) + ZT* N a 



d v g zz - 4,n z {x)s(y)5 e (z) + 2T a zz N a 



2n z {x)s{y) 



d v g zll - 2 nfl (x)s(y)5 e {z) - d»N z + 2T%N a 



N 



K ZVb 



(14) 



n^(x)s(y) 
N 



where K ab represents the extrinsic curvature except for the terms proportional to 5 e (z). Thus the gravitational 
equations on the codimension-1 brane at y = 0, off the codimcnsion-2 brane at z = 0, are given by 



1 2 M 4 

{5) G ab = —S ab + -^{Kab - Kg ab ). 



(15) 



Finally, let us derive the gravitational equations on the codimension-2 brane at y = z = 0. Plugging the 
components (fT4|) into equation (fT3]). and taking the integral over z across the codimension-2 brane, we obtain 
the localized terms to the codimcnsion-2 branc coming from the boundary contribution to the variations of the 
action Therefore, the boundary equations on the codimension-2 brane read 



n v {x) = 0, 



(16) 
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( %^^ + 2M (V -A; v ) + 4M r dz^n a (x){g"g^-g' IJi g' v )6 e {z). (17) 

4 4 4 J 

Equations ([TUf , (JTS]) and p7|) along with the bulk Einstein equations (J5J) are equations of motion in the cascading 
gravity model. It is obvious that when M| = 0, the DGP gravitational equations on the codimension-2 brane 
are recovered. 



3 Cosmology 
3.1 Bulk geometry 

In this section, we apply the formulation developed in the previous section to cosmology. We firstly present the 
bulk metric where a general Friedmann-Robertson- Walker (FRW) codimension-2 brane and a codimension-1 
brane are embedded. We start from the 6D Minkowski spacetime with the foliation 

ds 2 6 = (c 2 - fi 2 )dy 2 + (dz + (3dy) 2 -(l+{^-z + cy){H + ^j 2 dt 2 

+ a(tf(l + (-z + cy)H^5 ij dx i dx j , (18) 

where j3 and c are constants, a(t) is an arbitrary function of t and H{t) := —. "dot" denotes the derivative 
with respect to t. The metric (fT5)l clearly satisfies the 6D Einstein equation ([5]). Note that if j3 and c are time 
dependent the metric is not a solution to equation ([5]). On the y = z = hypersurfacc where the metric can be 
written as 

ds\ = - dt2 + a{tf5 l] dx l dx\ (19) 

the function a(t) can be interpreted as the cosmic scale factor of the FRW Universe, where H becomes the Hubble 
parameter. Imposing the ^-symmetry across the y — and z = hypcrsurfaces and choosing c 2 = 1 + /3 2 , the 
bulk metric reads 



dsl 



(l + /3 2 (1 - e{z) 2 ))dy 2 + (dz + Ps(y)e(z)dy) 2 (20) 
- (l + (4=2= + y/T+P\y\)(H + §)) 2 dt 2 + a 2 (l + (-^= + y/TTp\y\)H)X**dxi. 

From the comparison with the ADM metric in equation (TTJ), it is straightforward to read off the lapse function, 
shift vector and induced metric components. It is clear that the metric (|20[) is invariant under the parity 
transformation y — > —y or z — > —z. The assumption of regularized profile of e(z) is particularly important to 
see the contributions from the 6D bulk on the codimcnsion-2 brane. On the other hand, it is enough to assume 
the codimension 1-brane as a distributional object, where e(z) and 5 e (z) are treated as the usual sign function 
and the delta function, respectively. 

3.2 General cosmological equations on the codimension-2 brane 

We now derive the cosmological solutions in the 4D spacetime. Assuming metric (|2"0"]) , the induced metric on 
the y = hypersurfacc is given by 

dsj = dz 2 -(l+ ^ {Z)Z (H + -)) 2 dt 2 + a 2 (l+ ^ Z)Z h] 2 S i jdx i dx j . (21) 
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Calculating the extrinsic and intrinsic curvature tensors (see Appendix B for details) and using boundary 
equation (|17|) . we find the modified Friedmann equations on the codimension-2 brane 

MI 2 H 2 =p + PcS , (22) 
-M 2 (2H + 3H 2 ) =p + p cff , (23) 

where p and p are energy density and pressure of the matter localized to the codimension-2 brane, and p a g and 
p c H represent the effective energy density and pressure of the dark component 



(24) 



QMH f TT H\ , f , 05 J z) / 0e(z)z TT \ 2 

Poff := ; 5 [H+ — ) +4Mi / dz - H w 1+ -, = =H) 

The effective equation of state is given by 

Weil := (25) 



Poff 

i ^ 1 

H p c g 



^T+W J N /l + /3 2 (l-e(z)2) x /TT^V yrT^V 77 



Finally, the Bianchi identity on the codimcnsion-2 brane gives the nonconservation law 

a f P5Az) 0e(z)z f d ( H\/ 0e(z)z / fT> 

p + 3H(p + p) = AMi / ^ — H v ; , v ' -A 2—( + + JL AJ. — ( n + _ 

+ 3h(2 + -^^=(2H+^))\. (26) 

This may not be a surprising fact since the inclusion of a codimension-2 brane thickness is already somewhat 
beyond our theory ([2|) where the codimcnsion-2 brane is assumed to be an infinitcsimally thin object. In other 
words, to see the nontrivial bulk effect on the brane dynamics, it is necessary to include a finite brane thickness. 
For a de Sitter or Minkowski codimension-2 brane where H is constant the energy conservation law on the 
codimcnsion-2 brane is satisfied even including a finite thickness. 

From now on, in the main text of this paper we will focus on the cosmological dynamics on the codimcnsion-2 
brane. Concerning the dynamics on the codimension-I brane, we ask readers to sec the Appendix A. 

3.3 Small brane thickness approximation 

In general it is impossible to perform the integrals in equation (|2"41 analytically. We here use an approximation 
in which the codimension-2 brane thickness, tr, is much smaller than the size of the cosmic horizon H~ l . The 
details of the approximation are shown in Appendix C. We then find 



arctan(/3) + aC(0) ( H + — 



arctan(/3) + aC(0)H , (27) 



where we have defined 



C{0) := —jJ= + arctan(/3), (28) 

V 1 + 2 P 
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which is a non-negative even function of /3, vanishes at /3 = and monotonically increases toward £ as — > oo. 
The conservation law on the codimension-2 brane equation (|26p reduces to 

p + 3H(p + p)= Ul^C(P) ) + 4ffl . (29) 



Thus, due to a finite thickness of the codimension-2 brane there is always an energy exchange between the bulk 
and brane. The effective equation of state (|26|) becomes 



W eS = -1 - 772 SU^ I77T- — — • ( 3 °) 



2 7=f? ~ ^(arctan(/3) + f )) 



For \H\ < H 2 and since aH < f , 



2M 6 4 ctC(/3) 



_i _ g x/ 1 ^ 

^ 2 ^L-HM arctan(/3) - (31) 
4 De Sitter solutions and the features of degravitation 

In this section, we give the solutions with Minkowski and de Sitter codimension-2 branes in the cascading model. 
In particular, we discuss the possible features of degravitation. 

4.1 Minkowski codimension-2 brane 

The simplest solution is a Minkowski codimension-2 brane solution which is realized by setting H = and 
H = hi PU)| . In this case, the 6D bulk metric is given by 

(l + 2 (1 - e(z) 2 )y y 2 + (dz + f3s(y)e(z)dy) 2 - dt 2 + Sijdx l dx j . (32) 

The codimcnsion-2 brane is supported by the tension p = —p = A. Modified Friedmann equations (|22|) reduce 
to the single equation which relates the brane tension to the bulk geometry 

A = 4Mg arctan(/3), (33) 



ds 2 6 



s 



where the left-hand side shows the deficit angle in the bulk. For a Minkowski codimension-2 brane, K a b = 
and the codimension-1 brane geometry is also the 5D Minkowski spacetime. Since |arctan(/3)| < 1, the brane 
tension can be at most of order Mg. In Rcfs. [HIS], it has been shown that the cascading gravity is ghost free 
if the codimension-2 brane tension satisfies the following bound 

A > |^M 6 4 . (34) 

This condition is satisfied as long as two cross-over scales satisfy < 7-4 . Perturbations about the Minkowski 
codimcnsion-2 brane solution (|32[) and stability have been explicitly analyzed in the recent work Ref. [6]. 

4.2 De Sitter codimension-2 brane and degravitation 

Now we consider the de Sitter codimension-2 brane solution where a(t) = aoe Hot . In particular, we focus on 
the possible features of degravitation. In this case the 6D bulk metric (|20| reduces to 

ds 2 6 = (l + (3 2 {l- e(z) 2 ))dy 2 + {dz + (3s(y)e(z)dy) 2 - (l + ( ^1 + ^TT^\y\) H a f dt 2 
K 7 v ^l + p z J 

+ a 2 e 2Hat (l + {^^ + ^/l + ^\ y \)HX5 l3 dx l d^. (35) 

Here the codimension-2 brane can be supported only by the tension with p = —p = A. The effective cosmological 
equations (|22j) reduce to a single integral equation 



3M 2 H 2 = \ + PeS , (36) 



where 



Pes = -Poff = , H - 4M 6 / dz [ 1 + H j . (37) 

\/l + /3 2 J yJl + /P(l-e(z)z) V \/l + ^ 2 7 

In the following, we are going to obtain explicit de Sitter 3-brane solutions. The solution to equation (|36p 
gives the expansion rate of the de Sitter codimcnsion-2 brane. The solution of Hq > (< 0) represents an 
expanding (contracting) de Sitter Universe in terms of the flat slicing. In the case of A = 0, we discuss some 
self-accelerating solutions. 

4.2.1 Recovering DGP solutions 

In the absence of the bulk gravity where M@ — > 0, the modified cosmological equation on codimcnsion-2 brane 
reads 

3M 2 H 2 = X + -^M=Ho 7 (38) 

which corresponds to cosmological equation in DGP model with a codimension-2 brane tension. The solution 
is then given by 



" 3Mf 



#0 , 



where H^ > Hq '. If A > 0, irrespective of sign of (3 the solution of Hq + ^ represents the expanding de Sitter 
Universe and that of H Q represents the contracting de Sitter Universe. If — < A < 0, for /? > both 
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A 





a o 
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A > 


> 


Expanding 


Contracting 


< 


Expanding 


Contracting 


A = 


> 


Expanding 


Minkowski 


< 


Minkowski 


Contracting 




> 


Expanding 


Expanding 


/3 < 


Contracting 


Contracting 



Table 1: In this Table the classification of solutions is shown. The terms "expanding" and "contracting" denote 
the expanding and contracting de Sitter Universe, respectively. The term "Minkowski" denotes the Minkowski 
codimension-2 brane. Note that the same terms are used in the subsequent tables. 

the solutions of H^ represent the expanding de Sitter Universe while for B < both the solutions of i/g 
represent the contracting de Sitter Universe. If A = 0, for 8 > the solution of H}. ' = . 5 — gives the 

self-accelerating solution of the DGP model, while for (B < the solution of Hq = gives the normal branch 
Minkowski brane solution of DGP. But the self-accelerating solution in the DGP model is known to be unstable. 
The classification of solutions is shown in Table 1. 



4.2.2 Degravitation features in the limit of zero brane thickness 

In this section wc study the effective cosmological equations (|36[) in the limit of zero thickness of the codimension- 
2 brane. In this case equation ([56]) becomes 

fi M 3 B 

3717? Hi = \+ 51 H , A := A — 4Mi arctan(/3) . (40) 

y/l +/3 2 

The solution is simply given by equation (|39p classified in Table 1, with replacement of A with A. In the case of 
A = 0, for j3 > both the solutions of give rise to a self-accelerating Universe for Mf > 4M 4 M & (i+ff ) aictan p ^ 
For p < only the solution of fig provides the self-accelerating Universe, while leads to a contracting 

Universe. However, these self-accelerating solutions might be unstable against perturbations due to the possible 
existence of the ghost mode, because they do not satisfy the condition ([34]) . 

Now we argue the possible connections of our solutions with the degravitation expected in the cascading 
model. We focus on the case of A < and B > where both solutions represent the expanding de Sitter 
Universe. In the limiting case of 3M|/3 2 ^> —(1 + /3 2 )Af|A, wc have 



H (+) ~ 12/3 2 r 4 Mi 



« 1, (41) 



for -771 =0(1) and as long as r% <C r&. The physical meaning of equation (|4T1) is clearly understood as follows: 

2 8 

The Hq ' solution with the tension | — jjw 3 - < A < 4M| arctan /?, which is expected to be ghost free, gives a 
much smaller expansion rate than one in the self-accelerating branch in the 5D DGP model with an expansion 
rate of order . In addition, rewriting the Hq ^ solution in the same limit 

H ^'^jkO^<<^ xl (42) 

shows that the effect of the effective vacuum energy A into the Hubble expansion rate is suppressed by the ratio 
^ in comparison with the naive expectation from the ordinary 4D cosmology. This implies a deep connection of 
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our results with the dcgravitation idea. This result shows that cascading model can provide a mechanism that 
could support a small expansion rate in the presence of a large cosmological constant i.e., tension. Moreovere, 
although the suppression in the 6D model is not enough to explain the fine-tuning problem for obtaining a tiny 
expansion rate which is suggested from observations, the cascading gravity model may be extendable to higher 

M i_1 

dimensional cases, which may lead to more cross-over scales r, := , f *; +1 (i = 5, 6, 7, • • ■ ) with r% <C ri+i. Then 
we could expect that the fine-tuning problem is alleviated more than in the 6D model. 

4.2.3 Degravitation features with a small brane thickness 



In general, it is impossible to solve the integral equation ([36]) with (|37|) analytically when a finite thickness of 
the codimension-2 brane is taken into consideration. Here, we use the small thickness approximations given in 
the previous section and Appendix C. Using (|27[) . the equation (|36[) reduces to the quadratic equation for H$ 
which leads to the solutions 



H = H^ := 



3yi + £2 m| 



3MlP - 2 a/I + f3 2 aC(/3)M ( 



± \J (3Mf/3 - 2yjl + p 2 aC{p)Mlf + 3M|(1 + p 2 )\] . (43) 

There is no physical solution for 



- ~ _ {Mill - 2^TTW<?C{f3)Ml? (AA , 

X<K - 3(1 + WWi ' ( } 

The solutions are classified in Table 3. If A > the solution of Hq + ^ always represents the expanding de Sitter 
Universe. If A* < A < and /3 > 2 V 1+/3 ^^W M 6 both the solutions of represent the expanding de Sitter 
Universe. In the absence of the brane tension, A = 0, we can obtain self-accelerating solutions. If /3 > 0, for 
3M 5 3 ^ > 2^1 + /3 2 crC(/3)M 6 4 , both the solutions of H^ ±] give the self-accelerating universes, while if f3 < 
only the solution of Hq gives the self-accelerating universe. However, these self-accelerating solutions might 
be unstable against perturbations since they do not satisfy the condition (|34[) . 

Finally, we argue the connections of the Hq ' solution with the idea of degravitation. For j3 > and 
3M|/3 > 2^1 + f3 2 aC{/3)Mg, in the limit of A > A*, we obtain 

Hf 1 1 + I3 2 (n\ A 



and 



H (+) h _ ^2 12(3 2 \rj Ml < lf (45) 



3M| Vr 4 / ( _ 2^i+^aC{fi) ^ 



Hr'^^i-) : -M (46) 



3/3r 4 



for = 0(1) and ^ <C 1. This equation differs from equation (|4"2"j) by the factor (l — — ) 2 which 

depends on the thickness of the brane. Since r 4 3> <r, where r 4 could be macroscopic and a is microscopic, 
the effects of the brane thickness are small. Therefore, one could still conclude that the Hubble expansion rate 
is much smaller than the one in ordinary 4D cosmology, namely Hq ' 2 <C 3]qt|A|. Thus, irrespective of the 
inclusion of a small thickness, the cascading model exhibits features of degravitation. 
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P 


rr +) 




A > 


2v/l+/S 2 uCJ/ 6 4 

P > 3M.a 


Expanding 


Contracting 


3Mg 


Expanding 


Contracting 


A = 


_ 2\/l+(i2,CM 6 4 


Expanding 


Minkowski 


fl ^ 2Vl + /3 2 ^Ci\/| 
P < 3 A/ 3 


Minkowski 


Contracting 


A, < A < 


^ 2v/l+/S 2 »CJ/ 4 
P 3AJ| 


Expanding 


Expanding 


„ 2 V / 1 + /3 2 ^CAJ| 
" < ,ug 


Contracting 


Contracting 



Table 2: Classification of solutions with a codimension-2 brane thickness. 

5 Conclusion 

We have presented a formulation of the nonlinear effective gravitational theory in the 6D cascading gravity 
model. This model is a higher-dimensional extension of the 5D DGP braneworld model. In the simplest 6D 
model we are living on the codimcnsion-2 brane that is located on a codimension-1 brane embedded into a 
6D bulk. The particularly interesting expectations arc that this model may exhibit a degravitation where 
the gravitational force falls off sufficiently fast, and is also free from a ghost instability if the tension of the 
codimension-2 brane satisfies a bound. An important aim for presenting our formulation is to see whether in 
reality the idea of degravitation could work for the cascading model, through an explicit investigation of cosmo- 
logical solutions. The gravitational equations on the codimension-2 brane are composed of the contributions of 
the matter on the codimension-2 brane, the induced gravity on the codimcnsion-1 brane and the gravity in the 
6D bulk. The bulk contribution is given by integrating over the sixth direction across the codimension-2 brane 
along the codimcnsion-1 brane. After the derivation of the general equations of motion in the cascading model, 
we have applied them to cosmology where the codimcnsion-2 brane geometry is described by a flat Fricdmann- 
Robertson- Walker Universe, and obtained the modified Friedmann equations. In the zero thickness limit of the 
codimension-2 brane, the bulk contribution becomes an effective cosmological constant whose sign depends on 
the sign of j3. A finite thickness, however, leads to an energy exchange between the bulk and codimension-2 
brane, except for the cases where the codimension-2 brane geometry is exactly Minkowski or de Sitter. On the 
other hand, there is no energy exchange from or into the codimension-1 brane in any case (see Appendix A). 

Finally, we have discussed the Minkowski and de Sitter codimension-2 brane solutions. The Minkowski 
codimcnsion-2 brane solution is realized if both the bulk and codimcnsion-1 brane are empty and the codimension- 
2 brane tension takes a particular value determined by the bulk geometry. In the de Sitter brane solutions, the 
bulk gravity effect gives rise to a new branch of the solution which can give an expanding de Sitter codimension-2 
brane solution expected to be stable and has a much smaller expansion rate than in the original DGP model, 
which leads to alleviation of the fine-tuning problem. We have also shown some de Sitter solutions with the 
small expansion rate which could have deep connections with the idea of degravitation. This result shows that 
the cascading model provides a dynamical mechanism to resolve the cosmological constant problem. 
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A Dynamics on the codimension-1 brane 

Here we shall briefly mention the dynamics on the codimension-1 brane. The nonvanishing components of 
the extrinsic curvature tensors on the codimension-1 brane are given in the Appendix B. The components of 
the energy-momentum tensor on the codimension-1 brane S a b are determined by the gravitational equations 
([T3|) . From equation ([T3| and the fact of K zfJ _ = along the codimension-1 brane (// = t,i), we get S zfl ,=0. 
Thus, it is straightforward to check that the components of the extrinsic curvature tensor satisfy the condition 
Vb(A'^ — 8 b a K) = 0. Hence, the conservation law equation Vf,S 6 a = is also satisfied. Similarly in the case of a 
codimcnsion-2 brane without a thickness, since all the components of the 6D Einstein tensor vanish everywhere 
in the bulk, there is no energy exchange between the codimcnsion-2 brane and the bulk. But an inclusion of a 
finite thickness leads to an energy exchange between them if the codimension-2 brane geometry is not Minkowski 
or de Sitter symmetry. 

For a de Sitter codimension-2 brane, using equation (fT3|) together with (|56|) and ([57]) , the energy-momentum 
tensor of the matter localized to the codimcnsion-1 brane is given by 

(47) 



gz -6ffgA/ 5 3 + 8M 6 4 ff (yT+^ + pHpz) 

(^TT^ + f3H z) 2 

Thus the codimension-1 brane can be supported by the matter which has anisotropic pressure S z z ^ ^S* 1 ^. 
Note that they are regular in the codimcnsion-2 brane limit of z — >• 0. 



B Components of tensors 

B.l For a general cosmological brane 

On the codimension-1 brane, using metric (|20j) . we obtain the components of the extrinsic curvature K a b except 
for the contribution of the codimension-2 brane are given by 

1 s(y)(H + §) 



K % ■ 



N + {Pe{z)z + (1 + P 2 )\y\) (H + f ) 



N ^TTW+{Hz)z + {l + P 2 )\v\)H 
K zz - K zt = K zl = 0. (48) 

In the limit of y — > 0+, we obtain the nonvanishing components of the extrinsic curvature tensor except for the 
contribution of the codimension-2 brane 

K\ -> + " ^, K i j ^^=H S l , (49) 
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and hence the nonvanishing components of the combination appearing in the junction condition 



*W - -( E±i 3g Y (50) 

Given the 5D metric (|21[) . the nonvanishing components of the extrinsic curvature tensor on the codimension- 
1 brane (y = 0) are given by 

Vj = 6 i j (l + -?fik=H) - 13 He{z), K ti = 0. (51) 
In the codimension-2 brane limit of z — > 0+ 

K\ -> ; = + ><<', ' II. (52) 



The components of the five-dimensional Einstein tensor on the codimension-1 brane are given by 

12 _ 



^G z z = y 1)H ( 2H(H 2 + H)Az + 2H 2 + H 

(1 + HAz) (H + (H 2 + H)Az) 

(5) 3Q4 2 -l)ff 2 

' ~ (l + iMz) 2 ' 

(5) G\- = ?^ - ^ ; f3i?(i? 2 + ij)ylz + 2ij + 3iJ 2 V l 7, (53) 

(1 + HAz) (H + (H 2 + H)Az) V 1 

where we have defined A := . S 

On the codimension-2 brane, the nonvanishing components of Einstein tensor are given by 

(4) G* t = -3H 2 , Wffj = -S'^SH 2 + 2H). (54) 

The homogeneity and isotropy of the geometry on the codimension-2 brane restricts the form of the energy- 
momentum tensor to have only the diagonal components 

T\ = -p, V, /»>',. (55) 

B.2 For a de Sitter codimension-2 brane 

For the de Sitter codimension-2 brane the equations (f5T)]) reduce to 

K* v -V> v K = —=^° S^, K* Z -K = -^J^ , (56) 

where Hq is given by the solution of equation (|40j) . On the other hand, from equation (|53|) 

(5) G^ = 3g ° s6» v , (5) G% = 6g ° (57) 
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C Small thickness approximations 



Here we explain the small thickness approximation which has been used in the text. We use the representation 



e(z) = tanh I - I , which leads to S e (z) 



l,/ 
2 



e'(z) 



In the limits of a — > 0, they approach the usual 



2crcosh 2 ^ u 

sign s(z) and delta functions S(z), respectively. Thus a denotes the thickness of the codimension-2 brane. Here 
we assume crH -C 1, namely the brane thickness is much smaller than the size of the cosmological horizon, 
which is reasonable. Keeping the leading order corrections due to a finite thickness, for example the integral 
term in the effective dark energy density equation (|24j) becomes 



dz- 



^l + /? 2 (l-e(*) 2 
/3de 



1 



/3e(z)z 



v/l + Z^l-e 2 ) 



1 



2/3ae 2 



H 



(3de 



H 



i 

arctan(^) + [H 



1 



^(1- 
-§)<*?(/?), 



Pezje) 



H + — 



where C(/3) is defined in equation 
term in the effective pressure equation 



and shown in Fig. 1. The same procedure is also applied for the integral 



(58) 
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